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1, INTRODUCTION 


Program EAPROF (Emission Absorption Profiles) computes the 


| 
transverse profiles of infrared emission and extinction for an axisymmetric, 
axially uniform, cylindrical plume from radial profiles of gas temperature, 
pressure, and concentration and particle temperature and number density, 
The radiation model treats gas radiation transfer with band model methods 
and particle radiation transfer with the single-scattering approximation. 
The radiation model correctly couples the gas and particle components into 
a single emitting, absorbing and scattering medium, The program treats 
just one gas species and one particle species at a time. Development of 


the model is made in Ref, 1. 
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The gas band model is the Malkmus statistical model and employs 
either the Curtis-Godson (CG) or derivative (DR) approximations to 

: handle the inhomogeneity and nonisothermality of the plume. Lorentz, 
Doppler or Voigt line profiles may be used. 
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The single-scattering geometry used for particle radiation transport 
is shown in Fig. 1. The s-axis is the primary line of sight (LOS). (The LOS 
shown in Fig. 1 is the one that goes through the full plume diameter. As the | 
LOS is scanned out across the lateral extent of the plume, it cuts progressively 


shorter chords of the cylindrical plume), The c-axis is the scattering LOS, 


1, S. J. Young, Retrieval of Flow-Field Gas Temperature and Concentration 
in Low-Visibility Propellant Rocket Exhaust Plumes, U.S. Air Force 
Rocket Propulsion Laboratory, Edwards Air Force Base, Calif, (to be 
published). 
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Fig. 1. Single-Scattering Plume Geometry. 
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It is described by the value of s where it branches off the primary LOS 
and by the scattering angle 6 and the azimuthal angle ¢, The 
single-scattering approximation includes radiation emitted along the 
primary LOS and radiation that has been scattered once from the scattering 
LOS into the primary LOS. If the scattering LOS terminates on the nozzle 
exit plane, motor radiation scattered into the primary LOS is also included. 
The exit plane is modeled as a solid disc with uniform temperature and 
emissivity. 

Extinction of external radiation shown through the plume is assumed to 
be caused by gas absorption, particle absorption and particle outscattering. 
The single-scattering approximation does not allow inscattering to contribute 


to extinction calculations, 


The structure of the code is described in Section 2, and a brief 
description of the function of each subprogram is given, Preparation of 
input data for the code is described in Section 3. An example application 
for a plume containing H,0 and Al,03 as the gas and particle species, 
respectively, is given in Section 4. This example is taken from Ref, 1. 


A listing of the code is given in the Appendix. 
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Z. CODE STRUCTURE 


The organization of the code is shown in Fig. 2. EAPROF is the main 
program. The subroutine INPUT reads all data required for a run and 
processes it for compatability with the rest of the code. The function 
ZONEFIT interpolates on input radial profiles and fits the data on a grid 
of N equal thickness radial zones. The function ANGLFIT interpolates 
on the input differential scattering cross section and fits the data to the 


scattering-angle integration grid (also part of input). 


With the radial data fitted to a fixed radial grid, further interpolation 
is performed to fit the data to the primary and scattering lines of sight. 
These interpolations are performed by subroutines ZLOS and SLOS, 


respectively. 


QTHERM is the main routine for computing the thermal radiation 
source function along the primary or scattering LOS, This radiation arises 
from gas and particle emission. The function KDPARAM interpolates for 
gas band model parameters along a LOS from input tables of these data. 
PLANCK computes the blackbody function. The remaining functions listed 


under QTHERM are radiation functions that are variously employed depending 


on the lineshape and nonuniformity modes selected for the gas radiation band 


model, They are described in more detail in Refs. 2 and 3. 


2. S. J. Young, Description and Use of the Plume Radiation Code ATLES, 


SAMSO-TR-77-100, U.S. Air Force Space Division, El Segundo, Calif., 
13 May 1977. 


3. S. J. Young, Inversion of Plume Radiance and Absorptance Data for 
Temperature and Concentration, AFRPL-TR-78-60, U. S. Air Force 


Rocket Propulsion Laboratory, Edwards Air Force Base, Calif., 
29 September 1978, 
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Fig. 2. Program EAPROF and Associated Subprograms. 
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The single-scattering source function is computed in QSCAT by an 


integration over all scattering angles 6 and azimuthal angles ¢. Thermal 


radiation from the nozzle exit plane is also treated here for scattering lines 
of sight that terminate on this plane. 


The thermal source and scattering source functions are integrated over 
a LOS in subroutine RADNCE, The results are listed by OUTPUT. 
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The routines TRNSFR, INITIAL and STORE are bookkeeping routines. 


11 


—_——— — 


—-~ 


oo 


NE Ee SA Sl NS SO a gr 
emer g pt an sp ee eR eR A . —~ -— — ee ~ _ ea rei ae, regis _ J 


3, PREPARATION OF INPUT DATA 


A computational run of program EAPROF requires a set of program 
control cards to specify the mode of computation and to supply input data, 
Some program control cards simply specify a computation mode, some 
specify a computation mode and supply data, while others signal the code 
that blocks of auxiliary data are now to be read in. Each type of control 
card contains an alphanumeric name in the first ten card columns. These 
names must be spelled correctly and must be left-justified. If data are 
specified on a program control card, they must be entered in accordance 
with the format specification indicated in the detailed description of each 
card given below, All fields of the program control cards are 10 columns 
wide. In general, integer and alphanumeric data must be right-justified 
in their fields. Non-integer numerical data may be entered in either F or 


E formats (with decimal point and, for the latter, the exponential symbol E). 


_E-formatted data must be right-justified in their field. These same rules 


| apply to data entered on auxiliary card decks. The types of control cards 


and the data contained on them are illustrated in Fig. 3. A description of 


each type follows. 
ls Title Card - The card name is TITLE. Columns 11-80 of this 


card may be used for any identification title desired. 


Ze Calculation Data Card - The card name is CALCDATA,. 
SHAPE (format Al10) must be one of the alphanumeric values 
LORENTZ, DOPPLER, or VOIGT. INHOM (format A10) must 
have either the value CG (for the Curtis-Godson approximation) 
or DR(for the derivative approximation), NZONES (format I10) 
is the number of radial and transverse zones used by the spatial 
numerical integration routines. The maximum value of NZONES 
is 50. NSIGMA (format I10) is the number of equal-length segments 


that a scattering LOS is divided into for numerical integration, 
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If SFLAG (format J10) has the value 1, full gas-particle | 
calculations are carried out, If it has the value 0 (or blank), 


particle effects are ignored, 


Plume Data Card - The card name is PLMDATA, 


L1 (format E10) is the distance (cm) from the nozzle exit plane | 
to the observation scanning plane. L2 (format E10) is the | 
distance (cm) from the scanning plane to the end of the plume. 

TN (format E10) is the temperature (K) of the nozzle exit plane 


disc, and EN (format E10) is its emissivity, 


Species Card - The card name is SPECIES, The GDATA card 
described under card type 4 allows the entry of up to four different 
gas species, The variable GCOL (format I10) on the SPECIES card 
selects the desired species by assuming a value of 1 to 4, 
Similarly, PCOL (format 110) selects one of the three particle 
species read in by the control card type 8 by assuming a value 

of 1 to 3. 


Band Model Parameter Card - The card name is GPARAM. 
This card calls for the read-in of band model parameters for 
the gas species of interest. These parameters are listed if 
the variable PRINT (format Al10) has the value PRINT. The 


card deck structure for the parameters is given in Fig. 4. 


Gas Data Card - The card name is GDATA,. This card calls for 

the read-in of the radial profiles of gas temperature, pressure 

and concentration, The required deck structure is shown in Fig. 5. . 
If the variable PRINT (format Al0) has the value PRINT, these 


data are listed, 
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| 
7. Particle Parameters Card - The card name is PPARAM. | 


This card calls for the read-in of particle scattering data. 

The required deck structure is shown in Fig. 6. If the ; 
variable PRINT (format Al0) has the value PRINT, these data 
are listed. 


8. Particle Data Card - The card name is PDATA. This card calls 
for the read-in of the radial profiles of particle temperature 
and number density. The required deck structure is 
shown in Fig. 7. If the variable PRINT (format AI0) has the 
value PRINT, these data are listed. 


9. Angle Integration Grid - The card name is GRID. The card 
supplies and calls for the read in of data defining the grid over 
which angle integrations are carried out. NPHI (format I10) is 
the number of intervals over which the 360° azimuthal integration | 
is carried out. The scattering angle integration over 180° is 
carried out on the grid defined in Fig. 8. If the variable PRINT 
(format Al0) on the GRID card has the value PRINT, the 
scattering angle grid is listed. 


10, Execution Card - The card name is RUN, When this card is 
encountered, computations are begun using the data entered up 
to that point, and an output listing of the results is made. When 
the computation and results listing are completed, the program 
continues to read program control cards until a new RUN card 
is encountered. A new computation is then begun for all of the 
conditions and data of the first run except those which have been as 
changed by the intervening program control cards and auxiliary 
data decks. This process is repeated until an end-of-file card is 
encountered, With this feature, a large number of related runs 


can be made with one job submission. 
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Other than the requirement that all required data be specified before 
a RUN card is encountered and that auxiliary data immediately follow the 
control card that calls for them, the program control cards may be arranged 


in any order, 


Great care should be taken in the preparation of input data since very 
few checks of data consistency and setting of default values are provided, 
A general feature of data preparation is that, if particular data on a card 
are not required, they need not be specified. If none of the data on a control 


card is needed, that card need not be included. 
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4. EXAMPLE APPLICATION 


An example application of EAPROF is made here for a plume containing 
H,O as the active gas species and Al,O, as the particle species. <A discussion 
of this plume model is made in Ref, 1. The radial gas data are shown in Fig. 9. 
The particle radial profile is flat with the loading value N,, = 10°/em>, The 
he 
(Ref. 3) and are appropriate to a band center ~ 3985 cia and a bandpass of 


O band model parameters are the NERD wideband parameters of Fig. 10 


~ 300 cm”!, The nonresonant, self-broadening parameter is 
Yo = + 07394 cm! /atm. The efficiency for resonant self-broadening is 
6.53,and the efficiency for foreign gas broadening is 1. 00. 

Particle scattering cross sections were computed using Mie theory, the 
particle size distribution of Fig. 11 and the indexes of refraction 
m = 2.5] - 10.0018, 2.51 ~ 10.01 and 2.51 -i0.05. The first listed value is 
the accepted value for pure Al,O,. The stuff of real plumes is not likely to be 
so pure. The two other values are simply arbitrary values selected to 


parametrically study the problem. The value m = 2.51 - 10.0! is used in 


this example. The results for the scattering cross sections are shown in Fig. 12. 


Azimuthal integration was performed with a |6-point grid. The 
scattering angle integration grid and the manner in which it covers the 
weighting function sin 6 p(6) of the scattering source function integral 


[p(@) is a scaled value of do/dQ J] are illustrated in Fig. 13. 


Calculations were made for the Lorentz line shape and the CG 
approximation. The number of radial/transverse zones was 10, and scattering 
lines of sight were also divided into 10 segments for numerical integration. 
The distance from the exit plane to the observation plane was taken as 3 cm, 
and the distance from the observation plane to the end of the plume was fixed 
at 15 cm, The exit plane was modeled as a flat disc with uniform temperature 
T = 800 K and emissivity e = 0.75. 
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Radial Gas Data. 
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Fig. 12. Differential Scattering Cross Sections for Al,O,. 
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A listing of the input data for these calculation conditions is shown in 
Fig. 14. The resulting output is shown in Figs. 15 and 16, Figure 16 also 
shows the results for computations assuming gas-only and particle only 
conditions, The gas-only results were obtained by changing the SFLAG value 
on the CALCDATA card from | to 0. The particle-only results were obtained 
by setting SF LAG back to 1 and substituting a band model parameter card deck 
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in which the k values for all temperature were set to k = 10 cm7!/atm, 
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